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Influence of gap structures to speciflc heat in oriented magnetic flelds: 
Application to the orbital dependent superconductor, Sr2Ru04 

Hiroaki Kusunose* 

Department of Physics, Tohoku University, Sendai 980-8578 

We discuss influence of modulation of gap function and anisotropy of Fermi velocity to field 
angle dependences of upper critical field, Hc2, and specific heat, C, on the basis of the ap- 
proximate analytic solution in the quasiclassical formalism. Using 4-fold modulation of the gap 
function and the Fermi velocity in the single-band model, we demonstrate field and tempera- 
ture dependence of oscillatory amplitude of Hc2 and C. We apply the method to the effective 
two-band model to discuss the gap structure of Sr2Ru04, focusing on recent field angle-resolved 
experiments. It is shown that the gap structures with the intermediate magnitude of minima in 
[100] direction for 7 band, and tiny minima of gaps in [110] directions for a and f3 bands give 
consistent behaviors with experiments. The interplay of the above two gaps also explains the 
anomalous temperature dependence of in-plane Hc2 anisotropy, where the opposite contribution 
from the passive q/3 band is pronounced near Tc- 
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1. Introduction 

A determination of its gap symmetry has been a long 
standing issue in a course study of unconventional su- 
perconductivity. Because of a modulation of gap ampli- 
tude, thermodynamic properties at low temperature fol- 
low power-law behaviors, which give a hint of the struc- 
ture of the gap function.^' ^ However, this cannot de- 
termine an absolute direction of existing gap minima 
in principle. Moreover, it is rather difficult to distin- 
guish practically higher power-law behaviors due to point 
nodes from exponential behaviors of isotropic gap. Re- 
cently, a powerful method has been developed to mea- 
sure the modulation of gap amplitude directly. Namely, 
field angle dependences in oriented magnetic fields have 
shown oscillatory behaviors in thermal conductivity'^"^ 
and specific heat,^"^^ depending on detail structure of 
the gap function. 

In order to analyze the oscillatory behaviors, the so- 
called Doppler shift method has been used frequently, 
in which the energy shift in the low-energy quasiparti- 
cle (QP) spectrum due to supercurrent flowing around 
vortices is taken into account. ^^"-'^^ The method neglects 
the scattering by vortex cores and the overlap of the 



core states 
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It makes the theory to be valid in low- 
temperature and low- field limit, and it still remains un- 
known what is the range of its applicability in the (T, H) 
phase diagram. On the other hand, at fields near Hc2 
the analytic solution for spatially averaged quantities in 
the quasiclassical formalism is available with the rea- 
sonable approximation.^^' ^""^^ The present author ex- 
tended the approximate analytic method to anisotropic 
singlet and unitary triplet pairings as well as multiband 
superconductivity.^^ Comparisons with reliable numeri- 
cal calculations over all field range at the lowest tem- 
perature indicated that the method is competent for a 
semi-quantitative analysis of experiments in wide region 
of the (T, H) phase diagram. 



The purpose of this paper is to investigate the influ- 
ence of the gap structure and the anisotropy of the Fermi 
velocity to the oscillatory behaviors in Hc2 and C on the 
basis of the approximate analytic solution. With those 
knowledge we discuss the most plausible gap structures of 
the layered perovskite superconductor, Sr2Ru04, which 
has been surveyed extensively with use of the rotating 
magnetic fields.^-^O-^'^-^^ 

There now exists considerable evidence that Sr2Ru04 
pronounces the spin-triplet superconductivity (SC) with 
broken time reversal symmetry, for instance, the un- 
changed spin susceptibility in the Knight shift across 
Tc — 1.5 K,^^ the intrinsic magnetization detected by 
fiSR with the onset of SC,^^ and the absence of a coher- 
ence peak in l/TiT.^" 

Another important aspect of Sr2Ru04 is the orbital 
dependent superconductivity (ODS), which comes from 
three quasi two-dimensional cylindrical Fermi surfaces 
called a, (3 and 7. 31.32 r^j^^ latter is almost decoupled 
from the formers due to different parity of t2g orbitals 
with respect to the conducting layer. The ODS scenario 
can resolve the inconsistency between the power-law tem- 
perature dependences compatible with lines of tiny gap, 
and almost isotropic thermal conductivity under in-plane 
magnetic fields. Moreover, the ODS also explains the pe- 
culiar dependence of the penetration depth, which 
is incompatible with lines of tiny gap alone. ^"^ Recently, 
Deguchi and co-workers have measured the field-angle 
dependence of the specific heat to give a direct measure 
of thermally excited QP in accordance with the SC gap 
structure.^' They suggest that Sr2Ru04 exhibits the 
ODS and its primary SC gap in 7 band has minima in 
[100] direction, while the passive a, (3 bands have tiny 
gaps in [110] direction. 

Concerning the angle-resolved measurements, the 
anomalous temperature dependence of the in-plane Hc2 
anisotropy has been observed where its oscillatory ampli- 
tude with maximum in [110] direction {h]}^^^ = 1.5 T) 
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decreases at elevated temperature and even changes sign 
near Tc-^^ Although few theoretical arguments have been 
made for this puzzle in terms of the tetragonal symmetry 
breaking of the two-dimensional order parameter,^**' a 
direct evidence of the scenario has not been observed 
yet. It is important to recognize that a reasonable model 
should accompany correct field angle dependences with 
corresponding Hc2 anisotropy together with its anoma- 
lous temperature dependence. 

The paper is organized as follows. In the next sec- 
tion, we introduce our model of the gap function and 
the anisotropy of the Fermi velocity, then we give nec- 
essary formulas to discuss field angle dependence. In §3 
we demonstrate the influence of the gap structure and 
the anisotropy of the Fermi velocity to the field and the 
temperature dependence of oscillatory amplitude of Hc2 
and C in the single-band model, keeping the primary 7 
band of Sr2Ru04 in mind. In §4 we discuss the most 
plausible gap structures for Sr2Ru04 shedding a light on 
recent field angle-resolved experiments. The last section 
summarizes the paper. 

2. Model and formulation 

Let us begin with the model of the gap function. The 

d vector of the triplet pairing is factorized as d{k) = 
zA(p{k), where the angle dependence of the gap func- 
tion, ip{k) is normahzed as {\(p{k)\'^) = !,{■■■) being the 
average over the Fermi surface. Among various type of 
gap functions proposed so far we adopt the model pro- 
posed by Miyake and Narikiyo for 7 band,'^^ 



I sin^ [7ri? cos((/) — (jjo)] + sin^[7ri? sin((^ — 0o)] 



1 - Jo(27ri?) 

X sgn[sin(^- (/)o)], (2.1) 

where Jq{x) is the Bessel function of zcroth order, cf) is the 
azimuthal angle of fe, and (j)o denotes the position of one 
of the gap minima (we adopt (?io = 0, i.e., [100] direction 
for 7 band) . Although the original meaning of R is the ra- 
dius of the Fermi circle in unit oiir/a, a being the lattice 
constant, we regard it as a phenomenological parameter 
to characterize the magnitude of the gap minima. Note 
that the results in the present paper depends only on the 
magnitude of the gap, |</'(fe)|, so that there is no differ- 
ence between isotropic s wave and the chiral p-wave for 
example. The angle dependence of the gap magnitude is 
shown in Fig. 1. R = 1 corresponds to the gap with line of 
zeros, while i? = is the isotropic full gap. For compari- 
son, the dependence of / wave, <p(fc) = sin(20)e*'^, is 
also shown in Fig. 1. We will also apply the same form 
of eq. (2.1) to the secondary gap in the a/3 bands with 
the minima in [110] direction, i.e., (jjo = 7r/4. The mi- 
croscopic mechanism for these type of gap functions is 
discussed by Nomura and Yamada using the third-order 
perturbation theory in the 3-band Hubbard model. 

Next, we introduce the in-plane anisotropy of the 
Fermi velocity. The essential feature of the anisotropy 
with 4-fold symmetry is described by 




Fig. 1. The angle dependence of the gap magnitude for various 
R. The thin line represents the case of / wave for comparison. 




Fig. 2. The in-plane anisotropy of the Fermi velocity 



tV±X 



2 + 2a + a'^ 
{1 + a cos 2(j)) cos — a cos 2(j)) sin c 



■ (2.2) 



The anisotropy of v±{(j)) = y + is shown in Fig. 2. 

In the case of Sr2Ru04, a is small positive for 7 band, 
and almost isotropic, a ~ 0, for a/S bands. '^^ The mod- 
ulation of V along z direction is approximated roughly 
by Vz{k) = VzSgn{kz) for quasi two-dimensional Fermi 
surfaces. The two dimensionality is characterized by the 
ratio, X = {Vz)/{v'^) = '^v1/v\. The average of the Fermi 

velocity is defined as v — {v^ -\- Vy + v1) = \/ v']_ + . 

We now explain briefly the approximate analytic so- 
lution of the fundamental quasiclassical equations for 
weak-coupling superconductivity under magnetic fields. 
The solution can be obtained by the following approxi- 
mations: (i) the spatial dependence of the internal mag- 
netic field is averaged by _B, (ii) the vortex lattice struc- 
ture is expressed by the Abrikosov lattice, and (iii) the 
diagonal element of the Green's function is approximated 
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by the spatial average. Although all approximations men- 
tioned above is valid near iJc2, comparisons with reliable 
numerical calculations suggest that the solution is com- 
petent quantitatively in wide region of the (T, H) phase 
diagram except in very low T and H regions. For fur- 
ther detail we recommend the reader to refer the liter- 
ature.-^*' Here we only quote the expression of the 
free energy measured from the normal state, which is 
given for strongly type-II superconductors, S ~ ff, in 
the clean limit as. 



^sn/No = ^(A-i)«/ A^A^, + ^ne A| In ( 



T 



-A? In 



(2.3) 



where t is the band suffix (either 7 or a/3), A^o is the 
total density of states (DOS) in the normal state, ng = 
Nqi/Nq is the partial DOS of the band £, and a;„ = 
(2n -|- Vj-kT is the fermionic Matsubara frequency. Here 
A is the interaction matrix in imit of N^^ , 7 = 0.577216 
is the Euler's constant and the cut-off Wc/ttTc = 100 is 
used in the summation of the Matsubara frequency. The 
function is given by 



h 
with 



1 + 5/ 



2A 



t)^_L(fe) 



A2|<^,(fc)|2w^(m^„), (2.4) 



2A 



^l\^i{k)\^W\lU(,n), 



-1/2 



(2.5) 

where A = (2|e|If) is the magnetic length, uin = 
2Au}n/ve±{k), and W{z) = evic{—iz) is the Faddeeva 
function. Here Vi±{k) is the component of v perpendic- 
ular to the field, which is given for the in-plain field, 
H = H{cos(ph,sm^h,0) as 



ve±ik) 



2v^ 



2-|-2a- 



X ( sin((/) — (^/i) — acos20sin((^ -|- (/i/i) 



(2.6) 



Note that without the anisotropy in the velocity, a = 
0, the oscillatory behaviors depend only on the relative 
angle, (j>o — (j)h- In the single-band model the difference 
of X can be absorbed in the definition of Hc2 because v± 
appears only in the form oi A/v±. This is not the case for 
multiband models where Xi is generally different. Once 
we obtain the equilibrium free energy for given T and H 
by minimizing the free energy with respect to the gap 
magnitude A^, we can discuss the oscillatory behaviors 
of Hc2 and C. 

3. Results in single-band model 

Let us discuss the single-band model with the gap hav- 
ing minima in [100] direction {(po = 0). At T = the in- 
plane Hc2 anisotropy due to the modulation of the gap 
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Fig. 3. The in-plane field angle dependence of the upper critical 
field with line of zeros in [100] direction (o = 0, </>o = 0). 



and/or the anisotropy of v is given by the formula, 



H.2 = ^l^^ 



exp 



|^(fc)|2ln 




(3.7) 

From this expression the minima of Hc2 are obtained 
when the minima of the oscillation in the gap and in 
the perpendicular component of the velocity coincide. 
Thus, without the in-plane anisotropy in v, i.e., a = 0, 
the minima of Hc2 are simply realized when the field 
is parallel to the direction of the gap minima. On the 
other hand, without the gap modulation, i.e., i? = 0, the 
anisotropy of v gives rise to the minima of Hc2 when 
<j)h II [110] (<^h II [100]) for a > (a < 0). 

At finite temperatures we first discuss the case of the 
isotropic Fermi velocity, a = 0. Figure 3 shows the in- 
plane field angle dependence of Hc2 for i? = 1 at T/Tc = 
0, 0.1 and 0.2. As we expect it shows minima for (j)h, \\ (j)o 
(= and 7r/2). The oscillatory amplitude in what follows 
is defined as SA = A{(t)h \\ [110]) - A{(t)h \\ [100]) where 
A is either Hc2 or C. The temperature dependence of 
SHc2 scaled by H^2^^ ([100] denotes the field direction.) 
is shown in Fig. 4. The amplitude of the Hc2 oscillation is 
a monotonically decreasing function of T, and it becomes 
larger the smaller the size of the gap minimum is. The 
inset shows the temperature dependence of i?c2 for if || 
[100]. 

The field dependence of 6C at T/Tc = 0.1 is shown 
in Fig. 5. As the field decreases, the negative amplitude 
increases and changes sign at H*. It is emphasized that 
even for the gap with line of zeros, i? = 1, it reaches the 
maximum, and again decreases toward to zero. For larger 
modulation of the gap 6C has stronger field dependence. 
The overall behavior at the lowest temperature can be 
understood as follows. At higher fields the in- plane Hc2 
anisotropy yields the oscillation in H/ Hc2 since H is fixed 
in the experimental situation. It leads to the amplitude 
oscillation in A, which predominates the oscillatory be- 
havior in C. Therefore, the oscillatory amplitude of C 
has the opposite sign of 6Hc2- At lower fields H < H*, 
the Doppler shift argument can be applied. The lo- 
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Fig. 4. The temperature dependence of the oscillatory amplitude 
of Hc2 with different size of the gap minimum, i? (a = 0, 4n) = 0). 
The inset shows the temperature dependence of i?c2 for _ff || 
[100]. 



0.8 




Fig. 5. The field dependence of the oscillatory amplitude of the 
specific heat at T/Tc = 0.1 (a = 0, </)o =0). The inset shows the 
angle dependence of 6C. 

cal QP with the momentum k has the energy spectrum, 
Ef^ + Vs ■ k m the supercurrent flowing around vortices 
with the velocity Vs- This energy shift gives rise to a fi- 
nite DOS near the gap minima if |fc • Vs] > 1^(^)1- Since 
-L i?, the number of activated gap minima becomes 
largest in the case of ff || (po yielding a maximum in C . In 
the limit of ^ 0, 5C decreases again since the volume 
fraction of the locally activated QPs becomes small. 

It should be noted that the Doppler shift argument 
is valid only at low temperature. At elevated tempera- 
ture, positive contribution to 5C from the activated QPs 
near the gap minima becomes less pronounced and the 
negative contribution from the Hc2 anisotropy predom- 
inates. This is seen clearly by the field dependence of 
5C at T/Tc = 0.1, 0.2 and 0.3 as shown in Fig. 6. For 
T/Tc > 0.2 no sign changes occur in 6C. 

Next we discuss the effect of the anisotropy of v. For 
this purpose we consider the isotropic gap, R = 0. The 
temperature dependence of i5ifc2 with the anisotropy is 



Fig. 6. The oscillatory amplitude of specific heat at elevated tem- 
perature (a = 0, (/)o = 0). 
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Fig. 7. The temperature dependence of 5ifc2 for the isotropic full 
gap with the anisotropic Fermi velocity. 

shown in Fig. 7. The effect of anisotropy in (5iJc2 has 
the opposite sign of a and is monotonically decreasing 
(increasing) function of T for a < (a > 0). Note that 
the oscillatory amplitude due to the anisotropic velocity 
can be the same order of magnitude of that due to the 
gap modulation. 

Figure 8 shows the field dependence of 8C at T/Tc = 
0.1. Since there is no contribution from the Doppler- 
shifted QPs in this case (i? = 0), the oscillatory behaviors 
come purely from the anisotropy of -ffc2, which affects the 
oscillatory behaviors over weak field region. Thus, the os- 
cillatory amplitude of has the same sign of a and its 
magnitude decreases monotonically as H decreases. 

As will be discussed in the next section, the oscilla- 
tory behaviors in Sr2Ru04 is determined roughly by the 
contribution from the primary 7 band. Viewing a small 
positive anisotropy of v in the 7 band,"^* we expect that 
the primary gap in the 7 band has a modulation with 
an intermediate size of minima in [100] direction (i? ^ 
0.6-0.8, (/)o = 0) for consistency with experiments.^' 



J. Phys. Soc. Jpn. 



Full Paper 



Hirociki Kusunose 5 



T 1 T w. 1 1 1 1 r 




Fig. 8. The field dependence of SC for the isotropic full gap with 
the anisotropic Fermi velocity. 




Fig. 9. The temperature dependence of C at zero field. The open 

circle represents the experimental data taJjen from Ref.*^ The 
inset shows C(T) at H/h]}""'' = 0.1. 



4. Interplay of two-band superconductivity 

With the knowledge obtained in the previous section, 
let us discuss the oscillatory behaviors in Sr2Ru04. Since 
the a and /3 bands have similar characters of Fermi sur- 
face and the Cooper pair scattering between them is ex- 
pected not to be small, we adopt an effective two-band 
model for Sr2Ru04. We use the normal DOS, = 0.57 
and Hap = 0.43, based on the de Haas-van Alphen mea- 
surc;nic!nts.^^ The parameters of two dimensionality are 



estimated as x-y = 
respectively, and ( 



6.0 X 10"^ and Xa/3 = 3.0 x 10 



13 

4.0, from the experi- 



mental value of H^' /H)^' and the Fermi surface topol- 

ogy^** For the gap structure (pio) and the anisotropy 
of V (ai), we choose = 0.65, 



and a-v 



0.1 



for the primary 7 band. For the passive af3 bands, we 
adopt R.afj = 1, (t>apa = 7r/4 and ciafj = 0, which can ex- 
plain naturally the anomalous temperature dependence 
of 6Hc2 as will be shown shortly. 

Before discussing the oscillatory behaviors we first de- 
termine the interaction parameters for Sr2Ru04. This 
is done by the specific heat fitting ai H = 0. The best 
fitting is achieved by 



A 



Ai 
A 



A 

A2 



(4.8) 



with Ai = 0.4, A2/A1 = 0.65 and A/Ai = 0.12, as shown 
in Fig. 9. The open circle represents the experimental 
data taken from Ref."^^ The inset shows the results at 
H/h]^2°^ =0.1, which indicate that almost all the QPs 
in the a(3 band is already activated at this field. 

Figure 10 shows the temperature dependence of the 
in-plane Hc2 anisotropy. The behavior at low tempera- 
tures is similar to the single 7 band case in Fig. 4 ex- 
cept that the oscillatory amplitude is suppressed by the 
positive weak anisotropy in v. It is remarkable that the 
oscillatory amplitude changes its sign near Tc- This can 
be understood as an interplay of the primary gap in 7 
band having the minimum in [100] direction with the sec- 
ondary gap in Q!/3 band having the different minimum in 
[110] direction. The contribution from the secondary gap, 




Fig. 10. The field dependence of 5Hc2- The ratio of the two gaps 
at Hc2{T) is shown in the inset. 



which has the opposite sign of the primary gap's, is en- 
hanced near Tc yielding the negative amplitude in SHc2- 
The enhancement of the ratio of two gaps at Hc2{T) is 
shown in the inset of Fig. 10. This is the strong evi- 
dence that the secondary gaps in aP bands have the gap 
minima in [110] direction. The similar interplay of two- 
band superconductivity is observed in MgB2, in which 
the 3-dimcnsional passive tt band predominates over the 
2-dimensional primary a band showing the strong tem- 
perature dependence in H^^^ / H^^' 

Finally, we discuss the field dependence of SC for 
Sr2Ru04. The contributions from each bands are sep- 
arately shown in Fig. 11. Both the gap modulation with 
minima in [100] direction and the anisotropy in v with 
a > contributes to the positive oscillatory amplitude in 
SCj at low fields. Since the passive a/3 bands give no con- 
tribution to SC except in very low fields, total amplitude 
becomes smaller than the single-band result by a factor 
of n-y(=0.57). It is remarkable that the a0 contribution 
at low fields has positive sign. We expect naively that 
the gap modulation with minima in [110] direction gives 
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sign change near Tc- This is understood by the oppo- 
site contribution from the passive bands, which becomes 
predominant over the 7 band contribution near Tc- The 
observation that the interplay of two gaps is pronounced 
near Tc provides an alternative way to determine gap 
symmetry in passive bands, which is hardly accessible 
by ordinary experiments. With the assumed gap model 
the total behavior of 5C together with the a(3 band con- 
tributions is consistent with experiments. 

Acknovirledgment 

The author would like to thank T.M. Rice, M. Sigrist, 
I. Mazin, Y. Matsuda, K. Izawa, T. Watanabc, N. Nakai, 
K. Deguchi, T. Nomura, K. Kubo and Y. Kato for fruitful 
discussions. 



Fig. 11. The field dependence of SC. 



the negative sign of SCafi at low fields. However, this 
naive expectation based on the Dopplcr shift argument 
is valid only at very low temperature. Since the passive 
gap is small as compared with A^(T) even at low fields, 
T /Tc = 0.1 is regarded effectively in high temperature re- 
gion. Therefore, the positive contribution from a/J bands 
results as the case of T/T^ > 0.2 shown in Fig. 6 (note 
that ^0 = "■74 in this case). The overall behaviors of SC 
is consistent with experiments.^' ^° 

5. Summary 

We have investigated the influence of the modulation 
of the gap and the anisotropy of the Fermi velocity to 
the oscillatory behaviors oiHc2 and C, using the approx- 
imate analytic solution of the quasiclassical formalism. 

It is shown that for the single-band model the mini- 
mum direction of the in-plane iJc2 anisotropy coincides 
with either the minimum of the gap modulation or the 
Fermi velocity. If they are in different directions, the re- 
sultant amplitude becomes small due to the competing 
contributions from both effects. The absolute magnitude 
of the Hc2 anisotropy is a monotonically decreasing func- 
tion of T. Due to the Hc2 anisotropy the oscillatory be- 
havior of C at high fields shows the opposite tendency of 
the Hc2 anisotropy and its magnitude increases toward 
Hc2 ■ At low fields and low temperatures the effect of the 
so-called Doppler shift plays an important role yielding 
the opposite oscillatory behaviors to the high- field ones. 
Therefore we can safely conclude that the source of the 
oscillatory behaviors is the modulation of the gap, if a 
sign change is observed in the oscillatory amplitude of 
C (H) . It is noted however that the effect of the Doppler 
shift is pronounced only at sufficiently low temperatures. 

With a combination of the results above, we have dis- 
cussed the gap symmetry for Sr2Ru04. The oscillatory 
behavior of C is roughly determined by the primary gap 
in the 7 band. To be consistent with experiments we 
conclude that the primary gap has an intermediate mag- 
nitude of minima in [100] direction. Assuming that the 
secondary gap in a/3 bands has line of zeros in [110] di- 
rection we have shown that the in-plane Hc2 anisotropy 
exhibits the anomalous temperature dependence with a 
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